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We consider electron-phonon (e-ph) energy loss rate in 3D and 2D multi-component electron 
systems in semiconductors. We allow general asymmetry in the e-ph coupling constants (matrix 
elements), i.e., we allow that the coupling depends on the electron sub-system index. We derive a 
multi-component e-ph power loss formula, which takes into account the asymmetric coupling and 
links the total e-ph energy loss rate to the density response matrix of the total electron system. 
We write the density response matrix within mean field approximation, which leads to coexistence 
of symmetric energy loss rate Fs(T) and asymmetric energy loss rate Fa(T) with total energy loss 
rate F(T) = Fs(T) + Fa(T) at temperature T. The symmetric component Fs(T) is equivalent to 
the conventional single-sub-system energy loss rate in the literature, and in the Bloch-Griineisen 
limit we reproduce a set of well-known power laws Fs(T) oc T ns , where the prefactor and power 
ng depend on electron system dimensionality and electron mean free path. For Fa(T) we produce 
a new set of power laws Fa(T) oc T" a . Screening strongly reduces the symmetric coupling, but the 
asymmetric coupling is unscreened, provided that the inter-sub-system Coulomb interactions are 
strong. The lack of screening enhances Fa(T) and the total energy loss rate F(T). Especially, in 
the strong screening limit we find Fa(T) 3> Fs(T). A canonical example of strongly asymmetric 
e-ph matrix elements is the deformation potential coupling in many- valley semiconductors. 
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I. INTRODUCTION 

Electron-lattice energy loss in various bulk and low di- 
mensional semiconductor systems has attained a great 
deal of interest during the last few decades i One reason 
for the extensive studies of electron-phonon (e-ph) en- 
ergy relaxation arises directly from device applications. 
Another driving force is that the energy relaxation is con- 
nected to the versatile physics of microscopic electron- 
phonon interaction. Furthermore, the total e-ph en- 
ergy loss rate can be probed in hot electron experiments, 
which serve as a test for electron-lattice interaction. 

The e-ph energy loss is highly sensitive to various pa- 
rameters of the coupled e-ph system. The most impor- 
tant parameter is of course the nature of the electron- 
phonon interaction. In single-valley semiconductors the 
electron-phonon interaction can be typically described 
by deformation potential (DP) or piezoelectric coupling 
constants. 2 In quantum wells also other types of coupling 
due to heterointerface^ or electric fields induced quantum 
confinement may be important. In the above-mentioned 
electron lattice coupling mechanisms the dependence of 
the e-ph matrix elements on the electronic variables, such 
as momentum, can be typically ignored. The situation is 
different, e.g., in metals where momentum dependency 
must be included due to the high Fermi level. 5 How- 
ever, by simply setting the e-ph matrix elements to a 
constant is not the whole story even for DP coupling in 
semiconductors. In many-valley (MV) systems, where 
the conduction band minimum consist of several equiv- 
alent valleys, certain strain components lift the valley 
degeneracy, which directly shows that the e-ph DP cou- 



pling depends on the valley indices^^ It is well-known 
that this leads to an important role of the valley degree 
of freedom in the elastic properties and ultrasonic at- 
tenuation of doped MV semiconductors^^ However, it 
was not until recently that the effect of MV-DP coupling 
was investigated in context of the e-ph energy relaxation 
as welli^S This brings us to the topic of this work: we 
present a theoretical discussion on e-ph energy loss rate 
in carrier systems that consist of different components 
(sub-systems), which can have different e-ph coupling 
constants/matrix elements. Here we follow a terminol- 
ogy where the coupling is referred to as asymmetric if 
the e-ph matrix elements of the sub-systems are different 
and if they are the same the coupling is referred as sym- 
metric. A typical example of asymmetric e-ph matrix 
elements is the MV-DP coupling. 

This paper describes the total steady state e-ph en- 
ergy loss rate (or power loss) in semiconductors in the 
presence of symmetric and asymmetric e-ph coupling. 
We will assume that the lattice and charge carriers are 
in separate internal equilibrium. Within this so-called 
temperature mode l 11 ' 12 the total energy loss rate can 
be quite generally described by a symmetric power loss 
formulaiii^iiii 5 . 

P = F(T e ) - F{T ph ), (I) 

where T e and T p h are the electron and phonon tem- 
peratures, respectively. The exact form of the energy 
loss rate (function) F(T) requires microscopic derivation, 
but it can be phenomenologically linked to the e-ph en- 
ergy relaxation time T e - P h- This time scale is roughly 
equivalent to the thermal RC time constant defined by 
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Ge-phCe, where G e - ph = 5F{T)/ST is the macroscopic 
electron-phonon thermal conductivity and C e is the elec- 
tron heat capacity. At low temperatures T~_ ph oc Tf 
and C e oc T e , which leads to the well-known power law 
F(T) oc T n (n = p + 2). This temperature dependency 
has been experimentally verified in various electron and 
hole systems, which exhibit different prefactor and power 
n, by utilizing low temperature carrier heating techniques 
(see, e.g., Refs. [IlliSIliliiliSIillillllll and references 
there in). 

Single carrier component is equivalent to symmetric 
e-ph coupling, which is typically the unquestioned as- 
sumption made in microscopic derivations of Eq. {!]) and 
T e-ph- However, the existence of the asymmetric coupling 
can have a large effect, because it is typically unscreened 
due to inter-sub-system Coulomb interactions. This is 
the case even if the total electron system would seem 
to be at the strong screening limit, where the screening 
wave vector k is larger than thermal phonon wave vector 
qr = kBT/hv s (v s is the velocity of sound). In MV semi- 
conductors the dependence of the DP coupling on valley 
indices is a source of strong asymmetric coupling and, 
due to lack of screening, this enhances the e-ph energy 
relaxation. 9,10 In these systems the carrier sub-systems 
(valleys) have a full spatial overlap, but modern micro 
and nanofabrication techniques enable also realization of 
artificial bilayers where two 2D electron gases (2DEG) 
form spatially separated sub-systems. Now if asymmetric 
coupling exists in closely spaced bilayers the e-ph energy 
loss rate can be similarly enhanced as in MV systems. 
Here we will touch e-ph energy loss rate in both MV 
semiconductors and bilayer systems. 

This paper is organized as follows: In Section Til Al we 
derive the multi-component version of Eq. (fT]) by utiliz- 
ing Kogan's approac h 11 ' 12 , which connects F(T) to the 
electron density response function. We will allow that the 
e-ph matrix elements depend on the electron sub-system 
indices, which introduces the multi-component nature. 
We will assume that phonons cannot couple the electron 
sub-systems directly. For example, in many-valley sys- 
tems such direct coupling is relevant only for high en- 
ergy phononSfii 2 - Then in Sections III Bl and III CI we write 
the density response function (here a matrix describing 
the multi-component electron system) within mean field 
approximation, which eventually leads to coexistence of 
a symmetric energy loss rate Fs(T) and an asymmet- 
ric energy loss rate Fa(T) with the total energy loss 
rate F(T) = F A {T) + F S (T). In Sect. iHIAl we di scuss 
about the general aspects of Fa(T) and in Sect. IIII Bl 
derive several analytical results in the Bloch-Griineisen 
limit. For Fs(T) we reproduce a set of analytical results, 
which show well-known dependency on the parameters of 
the electron system, such as, dimensionality and electron 
mean free path l e . For Fa (T) we find a different (new) set 
of results and show that if the asymmetric and symmetric 
coupling have similar magnitude then Fa{T) 3> Fg(T) 
(provided that k exceeds qr)- In Section Till Bl we also 
view the role of Fa{T) in experiments and finally con- 



clude and summarize in Sect. II V 



II. THEORY 

A. Electron-phonon energy loss rate 

We assume that a phonon with momentum hq and en- 
ergy hw cannot directly couple the electron sub-systems, 
in which case the e-ph interaction is described via ma- 
trix elements Ai q j, where I = 1,2, ... ,L refers to elec- 
tron an sub-system and L is the total number of sub- 
systems. However, we will allow disorder induced elastic 
coupling of the sub-systems. The electron-phonon inter- 
action Hamiltonian is now given by 



(2) 



where b q (b q ) is the phonon annihilation (creation) oper- 
ator. The electron density operator p q i = ^ fc c\_ q ; c fc ; , 

where c kl (c k j) is the electron annihilation (creation) 
operator in sub-system I. Following Refs. [11IIT2I the 
e-ph interaction will be considered as a perturbation 
Hamiltonian that will cause transitions from initial state 
\i, {n q }i} with energy Si to final state \ f, {n q }f) with en- 
ergy £f. In state \e,{n q } p ) index e refers to electronic 
states and {n q } p is a set of phonon occupation numbers. 
The transition rate from initial to final state Wf% is given 
by the golden rule formula 

W fi = ^ |(/, {n q } f \ H e _ ph \i, {n q } t )\ 2 5(8* - £ f ). (3) 

By substituting Eq. @ into Eq. performing an 

ensemble average over the initial phonon and electronic 
states, and summing over final electronic states we obtain 
the phonon emission and absorption rates: 



W em (q) 



W ab (q) 



u 

x(N q + l)S(E hf - hw) (4a) 

Y^™ i \( f \52i Mq < lp - q ' 1 



xN q 6(E itf + fuS). 



(4b) 



Here Wi is the weighting factor for the electron many- 
body state, the energy difference Eij = Ei — Ef and 
N q = (n q ). We assume that the phonon system can 
be described by a thermal distribution N q = Nx ph (ui) = 

[exp(hu> / ksTph) — 1] 1 at temperature T p h and by a well 



defined phonon dispersion u> 



(we ignore phonon 



renormalization and hot phonon effects^). The total e- 
ph energy loss rate per d-dimensional electron volume V e 
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is given by the energy balance equation 

= _L - [e-« w /» flr - - (1 - e-^' kBT ")N q 
e q 

X MqjC^miq, w)M q , m , (5) 
Lrn 

where the latter equality is obtained by utilizing Wf = 
Wi ex.p[(Ei — Ef)/kBT e ] (T e is the electron temperature). 
The correlator Ci m (q,u) in Eq. JSJ) is defined by 

C l<m (q,u)) = 2Tih^w f (f\p q .i\i) (i\p qm \f) 

x5(E fti +Huj) (6) 

and it is connected to the density response ma- 
trix xi,m(q,w) through standard fluctuation dissipation 
relation 2 ^ 

(1 - e-^' kBT °)C ljm {q,w) = -2W e lm{xi, m {q,u)}. (7) 

Substituting Eq. into Eq. (J3J) leads to microscopic 
definition of F(T) in Eq. Q: 

P = F(T e ) - F(T ph ) (8a) 
F(T) = ^wN T (u)2e\M^{-x(q,u)}M q e q . (8b) 

For the sake of clarity we will denote matrices by a hat. 
Here the response matrix xi :Tn (q,u)) = xi,m(q,oj) de- 
pends on the properties of the electron system. There- 
fore, if the temperature dependency of %(q,w) is rel- 
evant then F(T) should be more generally written as 
F(T,T e ), where T e arises from the temperature depen- 
dency of x(<Zi w ) = X,(q)W,T e ). We have introduced a 
useful matrix notation where the coupling matrix M q is 
related to the e-ph matrix element M. q ,i through relation 

M q ,i= (Mghmi, (9) 

i=x,y,z 

where e» are components of the phonon polarization vec- 
tor e q (\e q \ = 1). If M q: i = Mg.m for all I, m or L = 1 Eq. 
@ reduces to power loss formula of Refs. [UI1Q1 (The 
equivalence with the expressions of Sergeev et. al— can 
be obtained by utilizing the mean field response function 
given in Rcf. 23 .) The elastic inter- valley scattering in- 
duced energy loss rate derived in Ref. [li] is also a special 
case of Eq. flgj). 

The obtained power loss formula applies to any type 
of e-ph coupling mechanism that does not depend on the 
electronic variables (momentum) in a single sub-system. 
However, here we will mainly limit our studies to defor- 
mation potential coupling to bulk acoustic phonons (in a 



volume Vph), in which case the coupling matrix is given 

by 

a7 '=vS& §§> (10) 

where p is the mass density of the crystal. We follow 
the notations of Ref. [l^: the L x 6 DP matrix 5 con- 
tains the deformation potential coupling constants and 
matrix S is the displacement-strain conversion matrix. 
The DP matrix depends on the properties of the elec- 
tron system. The form of S follows from the relation 
e Q/ 3 = \{du a /d[3 + dup/da) = iqu q (qpe a + q a ep)/2, 
which couples the six symmetric strain components e a p 
to components of displacement it. With S this relation 
is simplified to e a p — iqu q ^2 <S , ct /3, 7 e 7 . By definition S 
is a 6 x 3 matrix, which depends only on unit wave vector 
components q a — q a /q , i.e., it depends on the direction 
of phonon propagation. In our notation M q also con- 
tains the e-ph form-factors. However, when deriving the 
analytical results (Sect. IIII Bp we mainly set these form 
factors to unity, which is a reasonable approximation if 
we are far from the threshold q?t = 1 (t is the thickness 
of the electron system) . 

B. Energy loss rate with mean field density 
response 

Next, we assume that the electronic system can be 
described with the response of the non-interacting sys- 
tem Xo(Qj w ) under the external field plus the induced 
field of all electrons, i.e., we use the standard mean field 
approach [random phase approximation (RPA)] 24 . The 
RPA density response and dielectric function, generalized 
to a multi-component system, is given by 

x(g.w) = xo(q,w)e~ l (q,w), (Ha) 

e{q,u) = 1-V%{q,u), (lib) 

where matrix V contains the interaction potentials. The 
elements Vij = Vy = V d {q)F l d 3 (q), where V d (q) = 
e 2 /4eo (2/q) d 1 is the d-dimensional (d = 2,3) Fourier 
transform of the Coulomb interaction and F l d J (q) are the 
Coulomb form factors (eo is the background semicon- 
ductor dielectric constant). Here we will consider two 
cases: (1) the 2D/3D sub-systems are fully overlapping 
(many-valley/band, MV) or (2) two 2DEG layers sepa- 
rated by distance zq (bilayer). These cases are covered 
with form factor F l J {q) = [8ij + (1 — <5y) exp(— qz )] 3 ~ 
(if 1/q greatly exceeds the width of the individual layers 
in 2D). 

To study the energy loss rate of the electron system we 
must invert the dielectric matrix, which is a formidable 
task for arbitrary L, V and xo (<Zj w). For two sub-systems 
(L = 2) the inversion is tractable and well-known from 
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bilayer Coulomb drag effect (see e.g. Ref. 25) . Further- 
more, it is easy to show that if Vxo = A + (1 — 6)1 where 
matrix A follows symmetry Aij — A^j for all i, j and k, 
then e~ 1 (q,cu) = 



b' 1 ^1 + A/[b- Tr(A))\ . This form 
covers, e.g., the case where all interactions are similar and 
Xa{q,ui) is arbitrary: Vy = V(q), A = Vx (q,uj) and 
6 = 1. Here we will mainly concentrate on a transpar- 
ent model where all intra-sub-system and all inter-sub- 
system dynamics are similar, respectively. In this case 
we have {xohj = Xd^ij + XodQ- - <%)■ We further define 
Xo =Xd+(L- l)Xod and xi = -Lxod when xo(<?> w ) is 
given by 



xo(q,u) = xo(g,^)i + xi(q,^)Q A . 



(12) 



Here Qa = 1— Qs with {Qs}ij = 1/L and these matrices 
have a useful property QiQj = &yQj. We can either use 
the properties of Q,: to invert the dielectric function or 
note that now we also have Vxo = A + 61. We find 



£g 1 {q-,u)Q s +E A 1 {q,u)Q A , 



(13) 



where the scalar dielectric functions are Eg(q,u>) = 1 — 
[Vu + (L-i)Vu]xo(q,u) and £ A {q,u) = 1- (V n - 
V12) [xo(q, w) + xi{q,u)]. We can recognize that these 
both have a similar form to the standard RPA dielec- 
tric function. They are actually well-known from bi- 
layer physics and the poles of £g\(q,uj) are related to 
symmetric (optical) and asymmetric (acoustic) plasmons, 
respectively! 26 i 27 

Now by using Eqs. fIT]) -(fT3" ]l we can divide F(T) [Eq. 
into symmetric {Fg) and asymmetric (Fa) terms: 



F(T) = F S (T)+F A (T), 



(14a) 



M 2 



F t {T) = ^^JV r H2Im{-x«} , p ~ 2 (14b) 



(14c) 



where i = S, A , xs = Xo(q,u) and xa = Xo{q,^) + 
Xi(<Z,w). We have introduced yet another notation: 
the effective matrix elements Mf, which are defined by 
quadratic forms and screened with the scalar dielectric 
functions Si(q,uj). The matrix elements compactly de- 
scribe the contributions of symmetric (Mj) and asym- 
metric [Mj^ e-ph coupling. Note that from the prop- 
erties of Qi it follows that the effective matrix elements 
obey Mf > and, therefore, also F,(T) > 0. 

If the e-ph interaction is mediated through DP 
coupling [Eq. [[TO])] then quadratic form H? = 
e^S^E^Qg^A^Seq can be interpreted as a square of the 
deformation potential "constant". Note that H, is not 
a constant in general case: it depends on the phonon 
polarization (arises from e q ) and direction of propaga- 
tion (arises from S). If the DP coupling is mediated 
by one isotropic dilatational coupling constant S, which 
is the same for all sub-systems, then we simply have 
Zl=E*\e q -q/\q\\ 2 5 i<s . 



C. Response function of non-interacting electrons 

In order to perform quantitative analysis we must 
know the response function of the non-interacting sys- 
tem xo(<7,^)- In the pure limit ql e 3> 1 the off-diagonal 
elements of xo(q,uj) are zero and for the diagonal ele- 
ments we use the standard zero-Kelvin expressions : 2 ^ 29 



Xo(<?,^) = 



|| [H (a+)-H(a-) 
hp. 

k 



■S= + G(a+,a-[ 



3D 
2D 



(15a) 



u> , o 2hw . 

l + i K d (al ) 

qv F e F 



(15b) 



where a± = uj/qvp ± q/2kp, H(x) = 2x + (x 2 — 
l)ln[(a;-l)/(x + l)] [ ln(z) branch: |Imln(z)| < tt ] 
and G(a+,a_) = -(a\ - l) 1 / 2 + [20(|a_| - l)0(-a_) - 
l](aL — l) 1 / 2 . v — v(ef) is the (dimensionality depen- 
dent) density of states at Fermi level £f > fcgT and vf is 
the Fermi velocity. Eq. (|15b[) is small-q approximation, 
where K d {x) = (v/2) d - 2 (1 - x) (1 - x) {d ' 3)/2 . 

We include the effect of (static) disorder by introducing 
a phenomenological transport relaxation rate 7 = 1/t, 
which is connected to Drude mobility by /j, = er/m. 
We distinguish between intra-sub-system and inter-sub- 
system scattering rate which we denote by 70 and 71, 
respectively. The total scattering rate is given by 1/t = 
7 = 7o + (L — 1)71. Now the components of Eq. lfT2|) in 
the diffusive (hydrodynamic) limit ujt, ql e <C 1 are given 
by the particle number conserving many-band response 
functions of Kragler and Thomas^ (see also Ref. 



xo(q,u) = -v 



iD q 2 



oj + iD q 2 ' 



xi(q,u) = Xo0?>w) 



27 



iD q 2 (u + 27) + iD Q q 



(16a) 
:,(16b) 



where Dq 



2 v Fle is the diffusion constant and 
7 = -j7iis the total elastic inter- valley/band scattering 
rate. Note that xo(?: w ) has the familiar diffusion pole 
form 2 ^, while the inter-band term Xi(q^ w ) exhibits also a 
relaxation pole which is slightly shifted towards finite q. 
The absorptive part of xa — Xoiq^) + Xi(q> w ); which 
enters to Fa(T), is 



ImxA 



v (7 + D g 2 ) uj 
ft + D q 2 ) 2 +u< 



(17) 



This simply states that in the diffusive limit there are 
two competing time scales which determine the energy 
relaxation: inter- valley scattering time I/7 7 and diffusion 
time l/Dgq 2 over the length scale 1/q. 
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FIG. 1: Schematic cross-section of a two-component system 
and illustration of destruction of screening under asymmetri- 
cally coupling field. The electron gases have thickness t and 
they are separated by a distance zo. Note that zo > de- 
scribes a bilayer, where as zo — describes a many-valley 
system, which requires separation in reciprocal space. The 
solid curves describe the external (vi xt ) and effective (wf") 
scalar coupling potentials in the plane of electron gas / = 1,2. 
The ± signs describe the charge build up due to the external 
(and effective) field in the two electron gases. See Sect. IIII Al 
for further details. 



III. DISCUSSION AND ANALYTICAL RESULTS 
A. General aspects 

As Fi{T) > holds always the symmetric and asym- 
metric energy loss rates in Eq. (1 14[) provide a separate 
additive energy relaxation channels. The appearance 
of two clearly distinguishable energy relaxation terms 
(channels) F$ and Fa follows from our symmetric choice 
of the response Xo{Q,w) [Eq. (p~2|) ] . If we would relax 
the symmetry of Xo{q,u>) more complicated cross terms 
with variable signs would also appear. However, here 
we will concentrate on the symmetric response, which 
captures the essential physics. The sufficient condition 
for finite asymmetric energy loss rate Fa is that asym- 
metric coupling Ma differs from zero. The magnitude 
of Fa depends on diagonal and off-diagonal elements of 
Xo(q,v) through \ A = Xo + Xi = Xd - Xod- Therefore, 
Fa > even if there is no direct coupling between the 
sub-systems, i.e., even if Xod — (xi = 0). This follows 
from an internal (dynamic) image charge effect which is 
due to Coulomb interaction between the sub-systems. 

The internal image charge effect is schematically de- 
picted in Fig. [1] and it quantitatively follows from the 
properties of matrices Qi in the inverse dielectric function 
of Eq. (JT2J) : let us assume that L is even and that the ex- 
ternal field couples to the electron system through (scalar 
potential) v ext (q,ui) with the components {v ext }i = 

v ext = v ext,S + v ext,A(_ x y ( m p ig ft v ext,S = v ext,A^ 

Here vf xt describes the (phonon) field interaction poten- 
tial to sub-system I. Now the effective potential felt by 
the full electron system is v e ^ = 2~ 1 (q,oj)v ext where 
v f f = el s ' l (q,w)v ext > s +e^ 1 (q,u)v ext ' A . Note that full 
asymmetry of v ext is of course not required. It is enough 
that vf xt = i) / v^ xt — v + Sv for one fixed m ^ I and 
Qa picks this particular component. 



If the interaction between the sub-systems is strong [ 
(Vii — V12) is small] then e^g, w) ~ 1 and this quenches 
the screening of asymmetric coupling constants and en- 
hances the total electron coupling to the external field. 
Thus, the effect described in this paper is most im- 
portant in the strong screening limit where the ther- 
mal phonon wave length 2ir/qT exceeds the electron sys- 
tem screening length 2tt/k. In this case we expect that 
F S (T) < F A (T), if e A {q T ^ T ) ~ 1 (wr = k B T/H) and 
M\ ~ Mg. Note that the asymmetric coupling also 
plays a role in a medium screening systems qr ~ k when 
Fs(T) ~ Fa(T). On the other hand, the assumption 
M\ ~ M| may not always hold. It holds if asymmet- 
ric and symmetric coupling are directly mediated by DP 
interaction like in many-valley semiconductors or bilayer 
systems with different DP coupling constants. But the 
existence of non-zero M\ may be due to vibration of 
heterointerfaces^. or e-ph form-factors, which may give 
M| > M%. Then if £ S (q T ,w T ) > £a{qt,ut) still 
holds we may have Fa{T) ~ Fs(T) and both symmet- 
ric and asymmetric terms contribute to total energy flow 
rate. We will inspect few different cases where either 
F S {T) < F A (T) or F A (T) - F S (T) in detail below. 



B. Analytical studies of different systems 

Table Q] shows asymptotic low temperature expres- 
sions of Eq. (fT"4f deep below the Bloch-Gruneisen limit 
qr = 2kF for different systems. Note that Fs(T) is cal- 
culated assuming strong screening, but in many cases 
Fa(T) is simply the weak/zero screening limit of Fs(T) 
(see the discussion above). Before going into further de- 
tails we will briefly describe the general guidelines how 
the formulas in Table [I] are obtained and how they should 
be interpreted. First of all, the sum over q in Eq. (|14b[) is 
converted into integral and linear phonon dispersion re- 
lations u> = uj q = v s q are assumed. The sound velocity v s 
refers to both transversal and longitudinal modes of the 
crystal and, therefore, all expressions must be summed 
over these modes (if coupling to both type of modes ex- 
ists). Here we mainly concentrate on DP coupling and 
use compact notation S| ^ = e q S^'E) Q s,A^Se q for the 
DP coupling constants [as discussed below Eq. (|14p]. 
In the pure limit qrl e 3> 1 the formulas for Fs^aIT) 
are obtained by utilizing Eq. (|15b[) . In the diffusive 
limit qrh <C 1 they are obtained by utilizing Eqs. (|16[) . 
For 2D electron system vector q in refers to 

the parallel component q n = qs'm9, i.e., in 2D we set 
x(q,oj) — ► x{q n ,Lu) (9 is the angle between q and the 
normal of the 2D electrons). The dimensionality also 
affects the Coulomb interaction [as described below Eq. 
(jllap ] and is eventually seen in the screening wave vector 

k = 2 [e 2 Lv/Aeo\ 1,{d ^ 1] ■ Note that in the 2D diffusive 
limit the Fa with ln(<7r) factors are derived with isotropic 
DP "E? A = Sq. For arbitrary analytical expressions 
cannot be found i^L 
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TABLE I: Symmetric Fs(T) and asymmetric Fa(T) contributions of energy loss rate of Eqs. (|14|l . The pure and diffusive 
categories are denned by qrle 3> 1 and qrl e <C 1, respectively. The Limits column gives additional assumptions used in 
calculating Fa- zo — refers to many- valley systems and zq > to bilayers. Temperature dependency comes from thermal 
phonon wave vector qr ~ -jf^-- The length scale l a = l e vp/v s d, the DP constants E 2 = e^S^S^ QiESe q and the dielectric 
function ea = 1 + kzo/2. The bracket {• ■ ■ ) stands for average over solid angle. The Fa with ln(gT) factors are derived using 
isotropic = E 2 ,. Both Fs,a(T) oc gj are normalized with -^r^-^Bn-i, where B n -i = F(n)£(n) = J dxx n ~ x / [exp(a;) — 1]. 
See Section Till Bl for further details. 



Category 



Fs(T) 



Fa(T) 



Fa/Fs 



Limits 



3D pure 



3D diffusive 



3 /c2 \ 8 



Is (Hi) q% 



( IT ) 



ST 



— Qrle, qrh > 1 
2g T « e , (g T Z s )-l » 1 



Zq T h,k B T/hj< 1 



2D pure 4^ (sin 0E%) q 1 ? 



q T 
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2D diffusive 



f Hgg^ln(g T / s ) 



^f; 3 (sin 2 0Hi)g« 

5^H§g| In (sa/ s ?t) 
^(sin 2 0cos 2 0~|)gS («Z s O ?t) 2 



(k/ s ) 



^q T l,qrls > 1 
2g T Z, (g^)" 1 » 1 
q T l,k B T/lvy < 1 
(eAgrL) > 1 
(sAgr/sT 1 > 1 



> 



> x 



We have divided the results in Table U into four cate- 
gory: 3D pure, 2D pure, 3D diffusive and 2D diffusive. 
These are described with common Fg(T), respectively. 
The symmetric component Fs(T) is equivalent to the 
conventional single-sub-system energy loss rate in the 
literature. Indeed, if we assume coupling only to longi- 
tudinal phonons all symmetric energy loss rates F$ (T) 
agree with the expressions given in Refs. Ill flllfl4lfl5l and 
references there in. Note, however, that for 2D diffu- 
sive case we have to set 1/kI s — in order to find 
equal expression with the Fs(T) of Sergeev et alJ^. 
This small discrepancy follows from the approximation 
scheme of Ref. [Ta which ignores dynamical effects in 
es(q,uj q ). The importance of the dynamical effects in 
the diffusive xo(q,uj) are determined by the threshold 
lu/DoQ 2 = 1 or equivalently by qxl s — 1, which in- 
troduces the length scale l s = l e VF/v s d. However, the 
imaginary part of 2D RPA response is proportional to 

sin 2 9/ 1 + I 2 sin 2 9 ( qsm.9 + k) 2 and, therefore, in the 
strong screening limit it is actually the parameter nl s 
which is important and not qrh- This reasoning ob- 
viously applies also for 3D case. Thus, both 3D and 2D 
diffusive limit Fs(T) in Table U are valid if nl s , n/qx 3> 1. 

Let us next inspect the asymmetric part of F(T). We 



first focus on MV systems that arc fully overlapping 
(zq = 0) . First notice that in pure cases the ratio 
F A /F S = (n/q T ) m where m = 4(2) in 3D (2D). Thus, 
if screening is strong F A ^> Fg, as we would intuitively 
expect. In the diffusive limit there are two competing 
internal relaxation mechanisms in the electron system, 
elastic inter-valley scattering and diffusion, as already 
discussed below Eq. (|17p . When diffusion dominates over 
inter-valley scattering y/j/jqrle ^> 1 and in the opposite 
situation y/j/jqrle <C 1- In contrast to Fg(T) now the 
magnitude of qrl s strongly affects the results. Note that 
actually all F A in MV systems depend only very weakly 
on dimensionality and when inter- valley scattering domi- 
nates (y/j/^qrle, ksT/fvy <C 1) the 3D and 2D rates are 
similar, because then either screening or diffusion plays 
no role. Despite the finite mean free path and inter-valley 
scattering rate still Fa 3> Fg holds in Table U Recently, 
experiments on n + Si films were performed — in the range 
where the energy relaxation should be described by inter- 
valley scattering induced e-ph energy loss rate: the fourth 
and ninth Fa (from top) in Table U which have Fa(T) cx 
^(5^)4 cx T e . This special case for F A {T) is also 
derived in Ref. [l(| but by using a phenomenological 
approach which utilizes "classical" phonon attenuation 
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rate. The experimental F(T) for the n + Si samples coin- 
cide with the Table |T] and ultrasonic attenuation data of 
Ref. §22 

In 2D bilayer systems the finite layer separation Zq re- 
sults in EA(q,u) ^ 1. If z is small (qrzo <C 1) then 
£,4(9,^) ~ 1 + kzo/2 (= ea)- So far we have assumed 
unity e-ph form-factors, but now we partially relax this 
assumption and allow phase factors, which arise from fi- 
nite zqi in the e-ph coupling matrix. The finite layer 
separation introduces exp[±iq z zo/2] « 1 ±iq z zo/2 factor 
in the first (+) and the second (— ) line of the DP matrix 
S, which gives rise to asymmetric coupling. It is easy to 
show that now EJQa^ — \^Zq"E*QsB, which leads to 
the bilayer (zq > 0) Fa in Table |U For bilayers, where 
the asymmetric coupling is induced by phonon form fac- 
tors, the ratio Fa/Fs ~ 1 and depends on the magnitude 
of the parameter kzq- 

If we would relax the symmetry of xo(Q- w ) [Eq. (fl2|) ] 
then more complicated cross terms in addition of Fs,a{T) 
would also appear in F(T), as was already pointed out 
in Sect. IIII Al Especially in bilayer systems the study of 
such terms is an interesting and important problem of its 
own, but it will not be discussed in detail in this paper. 
Note, however, that if (Vu — V12) ~ [ea{q,^) ~ 1] 
and xo(<7, w ) is diagonal then Eq. ([M]) approximately 
holds with Imxs.A — Im {Tr[xa(q, u>)]} (provided that 
k > qx). In this case the 2D Fa with zq — in Ta- 
ble Q] qualitatively describe closely spaced bilayers with 
different DP coupling constants. One such system is 
the 2DEG bilayer realized in a double AlAs quantum 
wel l 33 ! 34 , where the two electron gases are from conduc- 
tion band valleys with different symmetry (depending on 
sample parameters) . 3 - 2DEG bilayer systems can also be 
tuned between two component and single-component sys- 
tems by external gates. Note, however, that the gate elec- 
trode can affect/contribute the e-ph relaxation processes 
if the gate-to-2DEG distance is small (< q^ 1 ), because 
then the gate is simply a one "component" of the total 
carrier system. This effect can be present in all gated 
samples at least at sufficiently low temperatures* 3 ^. 

An interesting special case "bilayer" is a single quan- 
tum well with two populated sub-bands with energies Eq 
and Ei . Phonons vibrate the heterointerfaces of the well 
(,i.e., change spatially the quantum well width) which is 
a source of e-ph coupling. 3 This coupling is defined by 
SE n ~ 2E n du z /dz and if E t - E = AE > k B T it is a 
source of asymmetric coupling with 'E/ A = q^e^AE. Now 
we find (utilizing Table [J), e -g-j for a pure system Fa <x 

AE 2 ( C ° S 2 W / ) q 5 T and F A /F S cx (n/q T ) 2 (AE/E) 2 (we 



have replaced Sg by a dilatational DP constant S). As 
typically AE/S <C 1 applies, the role of asymmetric cou- 
pling in the case of interface vibration is important only 
if screening is very strong. 

IV. SUMMARY AND CONCLUSIONS 

We have discussed on general aspects of the energy loss 
rate induced by symmetric and asymmetric e-ph cou- 
pling in 3D and 2D multi-component electron systems. 
We derived multi-component version of Kogan's power 
loss formula [Eq. ©], which takes into account the e-ph 
matrix element dependency on the electron sub-system 
indices and links the total e-ph energy loss rate to the 
density response function (matrix) of the electron sys- 
tem. We adopted standard mean field approximation to 
find the density response function. This led to coexis- 
tence of symmetric energy loss rate Fs{T) and asym- 
metric energy loss rate Fa(T) with total energy loss rate 
F{T) = F S {T) + F A {T) [Eq. HUD]. 

For Fg(T) we reproduced a set of well-known low tem- 
perature power laws Fs(T) oc T" s , where the prefac- 
tor and power ns depends, e.g., on electron system di- 
mensionality and electron mean free path l e . For F A (T) 
we derived a different set of power laws Fa(T) oc T nA . 
Screening strongly reduces the symmetric coupling and, 
therefore, also Fs(T). Whereas, the asymmetric coupling 
is typically unscreened, which enhances Fa(T) and the 
total energy loss rate F(T). This enhancement is large if 
the asymmetric and symmetric coupling constants have 
similar magnitude and screening is important. Under 
these assumptions Fa(T) 3> Fs(T), which we quantita- 
tively proved also for many special cases. 

In many-valley semiconductors the deformation poten- 
tial coupling constants depend on valley indices, which is 
a source of strong asymmetric e-ph coupling. Our find- 
ings agree with recent hot electron experiments on doped 
many- valley semiconductor (n + Si). Furthermore, the ef- 
fect described here should be present in various bilayer 
systems. 



Acknowledgments 

This work has been partially funded by the Academy 
of Finland (project 205467, CODE) and by the European 
Union (contract number 034236, SUBTLE). J. Ahopelto 
and J. M. Kivioja are acknowledged for useful discus- 
sions. 



* Electronic address: mika.prunnila@v tt.fi| 

1 B. K. Ridley, Rep. Prog. Phys. 54, 169 (1991). 

2 J. Singh, Physics of Semiconductors and Their Het- 
ero structures (McGraw-Hill, New York, 1993). 

3 F. T. Vasko and V. V. Mitin, Phys. Rev. B 52, 1500 (1995). 



4 B. A. Glavin, V. A. Kochelap, T. L. Linnik, and K. W. 
Kim, Phys. Rev. B 71, 081305 (2005). 

5 F. S. Khan and P. B. Allen, Phys. Rev. B 29, 3341 (1984). 

6 C. Herring and E. Vogt, Phys. Rev. 101, 944 (1956). 

7 R. W. Keyes, in Solid State Physics: Advances in Research 



8 



and Applications, edited by F. Seitz and D. Turnbull (Aca- 
demic Press, New York, 1967), vol. 20, pp. 37-90. 

8 M. Dutoit, Phys. Rev. B 3, 453 (1971). 

9 T. Sota and K. Suzuki, J. Phys. C: Solid State Phys. 15, 
6991 (1982). 

10 M. Prunnila, P. Kivinen, A. Savin, P. Torma, and 
J. Ahopelto, Phys. Rev. Lett. 95, 206602 (2005). 

11 S. M. Kogan, Fiz. Tverd. Tela 9, 2474 (1962), [Sov. Phys.- 
Solid State 4,1813-1819 (1963)]. 

12 S. Das Sarma, J. K. Jain, and R. Jalabert, Phys. Rev. B 
41, 3561 (1990). 

13 Y. Shinba, K. Nakamura, M. Fukuchi, and M. Sakata, J. 
Phys. Soc. Jpn. 51, 157 (1982). 

14 P. J. Price, J. Appl. Phys. 53, 6863 (1982). 

15 A. Sergeev, M. Y. Reizer, and V. Mitin, Phys. Rev. Lett. 
94, 136602 (2005). 

16 Y. Kawaguchi and S. Kawaji, Jpn. J. Appl. Phys. 21, L709 
(1982). 

17 K. Hirakawa and H. Sakaki, Appl. Phys. Lett. 49, 889 
(1986). 

18 R. Fletcher, V. M. Pudalov, Y. Feng, M. Tsaousidou, and 
P. N. Butcher, Phys. Rev. B 56, 12422 (1997). 

19 R. J. Zieve, D. E. Prober, and R. G. Wheeler, Phys. Rev. 
B 57, 2443 (1998). 

20 R. Fletcher, Y. Feng, C. T. Foxon, and J. J. Harris, Phys. 
Rev. B 61, 2028 (2000). 

21 X. P. A. Gao, G. S. Boebinger, J. A. P. Mills, A. P. 
Ramirez, L. N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 
94, 086402 (2005). 

22 R. Kubo, Rep. Prog. Phys. 29, 255 (1966). 

23 J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 (1986). 

24 D. Pines and P. Nozieres, The Theory of Quantum Liquids, 
vol. I (W. A. Benjamin, New York, 1966). 



25 A. G. Rojo, J. Phys. Condens. Matter 11, R31 (1999). 

26 S. Das Sarma and A. Madhukar, Phys. Rev. B 23, 805 
(1981). 

27 G. E. Santoro and G. F. Giuliani, Phys. Rev. B 37, 937 
(1988). 

28 J. Lindhard, Kgl. Danske Videnskab. Selskab Mat.-Fys. 
Medd. 28, 1 (1954). 

29 F. Stern, Phys. Rev. Lett. 18, 546 (1967). 

30 R. Kragler and H. Thomas, Z. Phys. B 39, 99 (1980). 

31 The ln(gr) terms in Table [T] result from integral 

hj(x,S 2 A ) = f °° dyy* [exp(y) - l]' 1 ( ^-T+Jjt ), 
where j = (j = 1) for zq — (zo > 0). 
Excellent approximation for isotropic H 2 ^ 



E and x 



given by (x, H 2 ,) 



iE 2 /~ dytfl [exp(y) - l]" 1 {ln(4j/x - 2) - j(2 - 1/yx)}, 
where j/o is defined from condition {■ ■ ■ } = 0. Eventually 
we find Iij (a^H 2 ,) ~ ctijSoBi ln(a;), where ctij is a 
constant of the order of unity. 

32 Note that r iv in Ref . Iiol is related to 7 through 7 = l/2n v . 

33 K. Vakili, Y. P. Shkolnikov, E. Tutuc, E. P. D. Poortere, 
and M. Shayegan, Phys. Rev. Lett. 92, 186404 (2004). 

34 M. Shayegan, E. D. Poortere, O. Gunawan, Y. Shkolnikov, 
E. Tutuc, and K. Vakili (2006), cond-mat/0606158. 

35 In AlAs F(T) also contains piezoelectric contributions, 
which are easy to calculate by adding the piezoelectric cou- 
pling to matrix M q . 

36 Note that the gate electron and 2DEG temperatures can 
be different, which can lead to photon cooling effects [D. R. 
Schmidt, R. J. Schoelkopf and A. N. Cleland, Phys. Rev. 
Lett. 93, 045901 (2004)]. 



